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Abstract 

We  discuss  spacecraft  Doppler  tracking  in  which  Doppler  data  recorded  on  the  ground  are 
linearly  combined  with  Doppler  measurements  made  on  board  a  spacecraft.  By  using  the  four-link 
radio  system  first  proposed  by  Vessot  and  LevineW,  we  derive  a  new  method  for  removing  from  the 
combined  data  the  frequency  fluctuations  due  to  the  Earth  troposphere,  ionosphere,  and  mechanical 
vibrations  of  the  antenna  on  the  ground.  Our  method  provides  also  a  way  for  reducing  by  several 
orders  of  magnitude,  at  selected  Fourier  components,  the  frequency  fluctuations  due  to  other  noise 
sources,  such  as  the  clock  on  board  the  spacecraft  or  the  antenna  and  buffeting  of  the  probe  by  non 
gravitational  forces  l^l.  In  this  respect  spacecraft  Doppler  tracking  can  be  regarded  as  a  xylophone 
defector. 

Estimates  of  the  sensitivities  achievable  by  this  xylophone  are  presented  for  two  tests  of  Einstein’ $ 
theory  of  relativity:  searches  for  gravitational  waves  and  measurements  of  the  gravitational  red 
shift. 

This  experimental  technique  could  be  extended  to  other  tests  of  the  theory  of  relativity,  and  to 
radio  science  experiments  that  rely  on  high-precision  Doppler  measurements. 


INTRODUCTION 

Spacecraft  Doppler  tracking  is  the  most  sensitive  technique  to  date  for  measuring  distances 
and  velocities  of  objects  in  the  solar  system,  leading  to  information  on  masses  and  higher  order 
moments  of  gravity  fields  of  planets,  their  satellites,  and  asteroids!3*  41.  Doppler  measurements 
have  also  been  utilized  to  search  for  gravitational  waves  in  the  millihertz  frequency  region!5*  *1, 
and  for  placing  upper  limits  on  amplitudes  of  signals  characterizing  relativistic  effects!7*  *1.  These 
Doppler  observations,  however,  suffer  from  noise  sources  that  can  be,  at  best,  partially  reduced 
or  calibrated  by  implementing  specialized  and  expensive  hardware.  The  fundamental  limitation 
is  imposed  by  the  frequency  fluctuations  inherent  in  the  docks  referencing  the  microwave 
system.  Current  generation  hydrogen  masers  achieve  their  best  performance  at  about  1000 
seconds  integration  time  with  a  fractional  frequency  stability  of  a  few  parts  in  10-16.  This 
integration  time  is  also  comparable  to  the  propagation  time  to  spacecraft  in  the  outer  solar 
system. 
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The  frequency  fluctuations  induced  by  the  intervening  media  have  severely  limited  the  sensitivities 
of  these  experiments.  Among  all  the  propagation  noise  sources,  the  troposphere  is  the  largest 
and  the  hardest  to  calibrate  to  a  reasonably  low  level.  Its  frequency  fluctuations  have  been 
estimated  to  be  as  large  as  10-13  at  1000  seconds  integration  timet9!. 

In  order  to  systematically  remove  the  frequency  fluctuations  due  to  the  troposphere  in  the 
Doppler  data,  it  was  pointed  out  by  Vessot  and  Levine!1!  and  Smarr  et  oU1#l  that  by  adding 
to  the  spacecraft  payload  a  highly  stable  frequency  standard,  a  Doppler  readout  system,  and 
by  utilizing  a  transponder  at  the  ground  antenna,  one  could  make  Doppler  one-way  (earth-to- 
spacecraft,  spacecraft-to-earth)  as  well  as  two-way  (spacecraft-earth-spacecraft,  earth-spacecraft- 
Earth)  measurements.  This  way  of  operation  makes  the  Doppler  link  totally  symmetric 
and  allows  the  complete  removal  of  the  frequency  fluctuations  due  to  the  earth  troposphere, 
ionosphere,  and  mechanical  vibrations  of  the  ground  antenna  by  properly  combining  the  Doppler 
data  recorded  on  the  ground  with  the  data  measured  on  the  spacecraft.  Their  proposed  scheme 
relied  on  the  possibility  of  flying  a  hydrogen  maser  on  a  dedicated  mission.  Although  current 
designs  of  hydrogen  masers  have  demanding  requirements  in  mass  and  power  consumption,  it 
seems  very  likely  that  by  the  beginning  of  the  next  century  new  space-qualified  atomic  clocks, 
with  frequency  stability  of  a  few  parts  in  10“ 16  at  1000  seconds  integration  time,  will  be  available. 
They  would  provide  a  sensitivity  gain  of  almost  a  factor  of  one  thousand  with  respect  to  the  best 
performance  crystal-driven  oscillators.  Although  this  clearly  would  imply  a  great  improvement 
in  the  technology  of  spacebome  clocks,  it  would  not  allow  us  to  reach  a  Doppler  sensitivity 
better  than  a  few  parts  10-16.  This  would  be  only  a  factor  of  five  or  ten  better  than  the  Doppler 
sensitivity  expected  to  be  achieved  on  the  future  Cassini  project,  a  NASA  mission  to  Saturn, 
which  will  take  advantage  of  a  high  radio  frequency  link  (32  GHz)  in  order  to  minimize  the 
plasma  noise,  and  will  use  a  specially  built  water  vapor  radiometer  for  calibrating  up  to  95% 
of  the  frequency  fluctuations  due  to  the  troposphere!11!. 

In  this  paper  we  adopt  the  radio  link  configuration  first  envisioned  by  Vessot  and  Levine  Ul, 
but  we  combine  the  Doppler  responses  measured  on  board  the  spacecraft  and  on  the  ground 
in  a  different  way,  as  it  will  be  shown  in  the  following  sections.  Furthermore  our  technique 
allows  us  to  reduce  by  several  orders  of  magnitude,  at  selected  Fourier  components,  the  noise 
due  to  the  clock  on  board  the  spacecraft,  or  to  the  antenna  and  buffeting  of  the  probe  by  non 
gravitational  forces!2!.  This  experimental  approach  could  also  be  extended  to  other  tests  of 
the  relativistic  theory  of  gravity  and  to  radio  science  experiments  that  rely  on  high-precision 
Doppler  measurements. 


DOPPLER  TRACKING  AS  A  XYLOPHONE  DETECTOR 

In  Doppler  tracking  experiments  a  distant  interplanetary  spacecraft  is  monitored  from  earth 
through  a  radio  link,  and  the  earth  and  the  spacecraft  act  as  test  particles.  In  a  one-way 
operation  a  radio  signal  of  nominal  frequency  Vq  referenced  to  an  on  board  clock  is  transmitted 
to  earth,  where  it  is  compared  to  a  signal  referenced  to  a  highly  stable  clock.  In  a  two-way 
operation  instead  a  radio  signal  of  frequency  vo  is  transmitted  to  the  spacecraft,  and  coherently 
transponded  back  to  earth.  In  both  configurations  relative  frequency  changes  Au/vo  as  functions 
of  time  are  measured,  and  the  physical  effects  the  experimenter  is  tiying  to  observe  appear  in 
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the  Doppler  observable  as  small  frequency  changes  of  well-defined  time  signature!3’  4»  5k 

If  a  Doppler  readout  system  is  added  to  the  spacecraft  radio  instrumentation,  and  a  transponder 
is  installed  at  the  ground  station  (Figure  1),  one-way  as  well  as  two-way  Doppler  data  can  also 
be  recorded  on  board  the  spacecraft!1!.  If  we  assume  the  earth  clock  and  the  on  board  clock 
to  be  synchronized,  then  the  one-way  and  two-way  Doppler  data  measured  at  time  t  on  the 
earth  (i?i(i),  E2{t)  respectively),  and  the  one-way  and  two-way  Doppler  measured  at  the  same 
time  t  on  the  spacecraft  (Si(t),  S2(t )),  have  the  following  analytic  expressions 


Ei  (t) 

= 

+  Cac(t  —  L)  —  CE(t)  -|-  T(t)  +  B(t  —  L)  + 

Asc(t  • 

-L)  +  RL&it)  +  P^t) 

(1) 

m) 

=  ^(t) 

+  CE(t~2L )  -  CE(t)  +  2 B(t~L)  + 

+T(t 

-2L)  +  T(t)  +  AE(t-2L)  +  Aac{t-L)  + 

+TRsc{t  —  L)  +  EL^t)  +  Pgj(t) 

(2) 

Si(t) 

=  s?W 

+  GE(t  —  L)  —  Cac(t)  -f-  T(t  —  L)  -f-  B(t)  -f- 

+A£(t  —  L)  +  ELsi(t)  +  PSl  (t) 

(3) 

S2(t) 

=  s$(t) 

-t~  Eac(t  —  2L)  —  Cac(t )  +  2T(t  —  fr)  + 

+B(t 

—  2L)  +  B(t )  +  Aac(t  —  2L)  +  AE(t  —  L)  + 

-\-TR£{t  —  L)  +  ELg^it)  +  Psj(t) 

(4) 

where  E%(t),  F^(t),  Sf(t),  S!J(t)  are  the  contributions  of  a  signal  (a  gravitational  wave  pulse 
or  a  red  shift  for  instance)  to  the  four  Doppler  data,  and  L  is  the  distance  to  the  spacecraft 
measured  in  seconds. 

In  Equations  (1,  2,  3,  4)  we  have  denoted  with  CE(t)  the  random  process  associated  with 
the  frequency  fluctuations  of  the  clock  on  the  earth,  Cac(t )  the  relative  frequency  fluctuations 
of  the  clock  on  board,  B{t)  the  joint  effect  of  the  noise  from  the  antenna  and  buffeting  of 
the  spacecraft  by  non  gravitational  forces,  T(t)  the  joint  frequency  fluctuations  due  to  the 
troposphere,  ionosphere,  and  ground  antenna,  AE(t)  the  noise  of  the  radio  transmitter  on  the 
ground,  Asc(t)  the  noise  of  the  radio  transmitter  on  board,  TRac(t)  and  TRE(t)  the  noise  due 
to  the  transponder  on  board  and  on  the  ground  respectively,  EL^(t),  ELe^  (t),  ELg^t),  and 
ELsiit)  the  noises  from  the  electronics  at  the  ground  station  and  on  the  spacecraft  in  the 
one-way  and  two-way  data,  and  Pe,  (t),  Pg^t),  Psl  (t),  and  the  frequency  fluctuations  due 
to  the  interplanetary  plasma.  Note  that  the  noise  due  to  the  transmitters  on  the  ground  and 
on  board  have  been  denoted  with  the  same  random  processes  (A£(t)  and  Aac(t)  respectively) 
in  the  four  Doppler  responses.  This  is  correct  as  long  as  the  two  radio  signals  of  frequencies 
vo  and  Vo  transmitted  from  the  earth  and  the  spacecraft  are  amplified  within  the  operational 
bandwidth  (typically  forty  to  fifty  megahertz)  of  the  same  transmitters!13’  14l.  The  Doppler  data 
Si(t )  and  S2(t)  are  then  time-tagged,  and  telemetered  back  to  earth  in  real  time  or  at  a  later 
time  during  the  mission. 

In  Equations  (1,  2,  3,  4)  it  is  important  to  note  the  characteristic  time  signatures  of  the 
noises  CE(t),  Cec(t),  B(t),  T(t),  Ag(t),  and  Aac(t)l9<  19»121.  It  was  first  pointed  out  by  Vessot 
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and  Levine  I1!  that,  by  properly  combining  some  of  the  four  Doppler  data,  it  was  possible  to 
calibrate  the  frequency  fluctuations  of  the  troposphere,  ionosphere,  and  ground  antenna  noise, 
T(t).  Their  pioneering  work,  however,  left  open  the  question  on  whether  their  existed  some 
other,  perhaps  more  complicated,  linear  combinations  of  the  data  that  would  further  improve 
the  sensitivity  of  Doppler  tracking.  In  what  follows  we  answer  this  question,  and  derive  a 
method  that  allows  us  to  uniquely  identify  an  optimal  way  of  combining  the  data. 

Let  Ei(f)  be  the  Fourier  transform  of  the  time  series  E\{t) 

Ei(f)  =  /  E.it)  e2**  dt  ,  (5) 

J— OO 

and  similarly  let  us  denote  by  E2(f),  Si(/),  S2(f)  the  Fourier  transforms  of  and 

S2(t)  respectively.  The  most  general  linear  combination  of  the  four  Doppler  data  given  in  Eqs. 
(1,  2,  3,  4),  can  be  written  in  the  Fourier  domain  as  follows: 

y(/)  =  a(/,L)  Ex{f)  +  b(/,L)  E2{f)  +  c(/,L)  Si{f)  Hr  d(/,L)  S2[f)  (6) 

where  the  coefficients  a,  b,  c,  d  are  for  the  moment  arbitrary  functions  of  /  and  L.  If  we 
substitute  in  Equation  (6)  the  Fourier  transforms  of  Eqs.  (1,  2,  3,  4)  we  deduce  the  following 
expression 

Vif)  =  [a  m)  +  b  E°(f)  +  c  5J(/)  +  d  Sf(/)]  + 

+Ci(f)  [-  a  +  b  (e4**^  -  1)  +  c  e2r'*L]  + 

+CTc(/)  [o  e2”*1  -  c  +  d  {e4*ifL  - 1)]  + 

+f(/)  [a  +  b  (e4iri/i  + 1)  +  c  e2rifL  +  2d  e2"^]  + 

+B(/)  [a  e2*ifL  +  2b  e2wifL  +  c  +  d  (eirifL  + 1)]  + 

+TbU)  [ft  e^L  +  c  e2”*1  +  d  e2^L]  + 

+A^c(f)  [o  e^L  +  b  e2^L  +  d  eMfL\  + 

+a[ETElif)  +  +  b[T%c(f)  e2"'L  +  ET^if)  +  /£(/)]  + 

+c[ELSl(f)  +  J£(/)]  +  d[fRE{f)e2^L  +  EL^f)  +  *£(/)]  .  (7) 

The  four  coefficients  a,  b,  c,  d,  can  be  determined  by  requiring  the  transfer  functions  of  the 
random  processes  CE(f),  Csc(f),  f(/),  B(/),  AE(f),  ^*c(/)  in  Equation  (7)  to  be  simultaneously 
equal  to  zero,  and  by  further  checking  that  this  solution  gives  a  non-zero  signal  in  the 
corresponding  combined  data.  This  condition  implies  that  a,  b,  c,  d  must  satisfy  a  homogeneous 
linear  system  of  six  equations  in  four  unknowns.  We  calculated  the  rank  of  the  (6  x  4)  matrix 
associated  with  this  linear  system  by  using  the  algebraic  computer  language  Matkematica,  and 
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we  found  it  to  be  equal  to  two.  The  corresponding  solution  can  be  written  in  the  following 
way 

«(/,  L)  =  c(f,L)  e~2*ifL  -  d(f,L)  [ewfL  -  e^L\ 

b(f,  L)  =  -  c(f,L)  e~2irifL  -  d(f,L)  e~2wifL  ,  (8) 

where  c  and  d  can  be  any  arbitrary  complex  functions  not  simultaneously  equal  to  zero.  After 
substituting  Equation  (8)  into  Equation  (7),  we  have  derived  the  expressions  for  the  signal 
in  the  case  of  a  gravitational  wave  pulse  or  in  the  variation  of  the  gravitational  potential  as 
experienced  by  a  spacecraft  orbiting  a  celestial  body  (redshift  measurements).  We  found  that 
for  both  these  signals  their  combined  Doppler  responses  (Equation  (7))  also  vanish.  These 
results  imply  that,  at  any  Fourier  frequency  /  and  for  these  two  specific  Doppler  experiments, 
we  can  remove  only  one  of  the  considered  noise  sources.  Among  all  the  noise  sources  affecting 
spacecraft  Doppler  tracking,  the  frequency  fluctuations  due  to  the  troposphere,  ionosphere,  and 
mechanical  vibrations  of  the  ground  antenna,  T(/),  are  the  largest.  We  can  choose  a,  b,  c,  d 
in  such  a  way  that  the  transfer  function  of  T(/)  in  the  combined  data  is  equal  to  zero.  From 
Equation  (7)  we  find  that  a,  b,  c,  d,  must  satisfy  the  following  equation 

a(f,L)  =  -  b(f,L)  [e^L  +  l]  -  c(f,L)  e™fL  -  2 d(f,L)  e2^L  .  (9) 

Since  6,  c,  d  can  not  be  equal  to  zero  simultaneously,  we  will  choose  c  to  be  equal  to  1/2,  and 
b,  d  to  be  equal  to  zero.  In  other  words  we  will  consider  only  linear  combinations  of  one-way 
Doppler  data.  Note  that  with  this  choice  we  eliminate  from  y(t)  the  frequency  fluctuations  due 
to  the  transponders  and  the  interplanetary  plasma  affecting  the  two-way  Doppler  data.  These 
considerations  imply  the  following  expression  for  y(/) 

5(f)  =  \  [«?(/)  -  e&f)'?’,1L]  +  Ci(J)  eW  -  i c£(/)  [«*"'* +l]  + 

+jB(/)  [l  -  e**4]  +  —  [Ae(/1  -  A„(/)  e2"«]  + 

[P*(J)  ~  Pk(f)  eM/£]  +  j  -  ELe,(/)  e2*1^]  .  (10) 

Equation  (10)  shows  that  the  transfer  functions  of  the  noise  of  the  on  board  clock,  CM), 
and  of  buffeting  B(f),  can  in  principle  be  set  to  zero  (not  simultaneously)  at  specific  Fourier 
frequencies.  In  searches  for  gravitational  wave  pulses  it  has  been  shown!2!  that  one  can  reduce 
by  several  orders  of  magnitudes  the  noise  of  the  on-board  clock  at  the  nulls  of  its  transfer 
function  without  removing  the  gravitational  wave  signal.  For  redshift  experiments  instead,  a 
cancellation  of  the  noise  of  the  on-board  clock  at  those  frequencies  removes  also  the  signal. 
If,  however,  measurements  of  y(t)  axe  made  at  the  Fourier  frequencies  for  which  the  transfer 
function  of  £?(/)  is  equal  to  zero,  some  further  improvement  in  sensitivity  can  be  achieved!15!. 
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EXPECTED  XYLOPHONE  SENSITIVITIES 


In  what  follows  we  provide  the  expression  for  the  noise  in  y(t)  at  the  xylophone  frequencies 
in  the  case  of  gravitational  wave  searches.  The  analogous  expression  for  redshift  experiments 
can  easily  be  derived  in  similar  fashion. 

Let  6  be  the  time  interval  over  which  a  Doppler  tracking  search  for  gravitational  waves 
is  performed.  The  corresponding  frequency  resolution  A /  of  the  data  is  equal  to  1/6.  This 
implies  that  the  fluctuations  of  the  clock  on  board  can  be  minimized  at  the  following  frequencies 


fk  = 


(2fc-l)  A/ 
4L  2 


J 


k  —  1)2, 3, .... 


(11) 


At  these  frequencies,  and  to  first  order  in  A /  L,  the  Doppler  response  y(fk)  is  equal  to 


am  «  \  [s?(/»)  +  i  (-i)fc  Eim]  +  <  (-i)*+1  cEm  ±  («a/l)  c„m  + 

+B(fk)  +  \  [ETSim  -  iEL+CM  M)*+1]  + 

+\  p£(A)  -  i  PEiUk)  (-l)t+ 

+\  [aTc m  +  i  Mfk)  (-i)fc+1' 


(12) 


For  a  typical  gravitational  wave  experiment,  8  —  40  days,  and  A /  =  3.0  x  10-7.  Therefore,  the 
frequency  fluctuations  of  a  clock  on  board  a  spacecraft  that  is  out  to  1  AU  are  reduced  at  the 
xylophone  frequencies  by  the  following  amount: 


ttA  fL 
c 


=  4.7 


x  10"4  . 


We  should  point  out,  however,  that  these  resonant  frequencies  in  general  will  not  be  constant, 
since  the  distance  to  the  spacecraft  will  change  over  a  time  interval  of  forty  days.  As  an 
example,  however,  let  us  assume  again  L  =  1  AU,  6  ~  40  days,  and  /  =  5  x  10-4  Hz.  The 
variation  in  spacecraft  distance  corresponding  to  a  frequency  change  equal  to  the  resolution  bin 
width  (3  x  10“7  Hz)  is  equal  to  1.0  x  105  Km.  Tfrijectory  configurations  fulfilling  a  requirement 
compatible  to  the  one  just  derived  have  been  observed  during  past  spacecraft  missions!9),  and 
therefore  we  do  not  expect  this  to  be  a  limiting  factor. 

From  Equation  (12)  we  can  estimate  the  expected  root-mean-squared  (r.m.s.)  noise  level  er(/t) 
of  the  frequency  fluctuations  in  the  bins  of  width  A/,  around  the  frequencies  fk  (fc  =  1,2,3, ....). 
This  is  given  by  the  following  expression 


°Uk)  =  [Sym  Af}1'2  ,  k  =  1,2,3, ....  , 


(13) 
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where  5y(/*)  is  the  one-sided  power  spectral  density  of  the  noise  sources  in  the  Doppler  response 
y(t)  at  the  frequency  /*.  In  what  follows  we  will  assume  that  the  random  processes  representing 
these  noises  are  uncorrelated  with  each  other,  and  their  one-sided  power  spectral  densities  are 
as  given  in  Table  I.  In  this  table  we  have  assumed  a  frequency  stability  of  1.0  x  10~16  at  1000 
seconds  integration  time  for  the  clock  at  the  ground  station.  Although  this  is  a  factor  of  four 
better  than  what  has  been  measured  so  fart11!,  it  seems  very  likely  that  by  the  beginning  of 
next  century  such  a  sensitivity  can  be  achieved.  As  far  as  the  other  sensitivity  figures  provided 
in  Table  I  are  concerned,  they  are  as  given  in  the  Riley  et  al.  report!11*.  This  document  is 
a  summary  of  a  detailed  study,  performed  jointly  by  scientists  and  engineers  of  NASA’s  Jet 
Propulsion  Laboratory  and  the  Italian  Space  Agency  (ASI)  Alenia  Spazio,  for  assessing  the 
magnitude  and  spectral  characteristics  of  the  noise  sources  that  will  determine  the  Doppler 
sensitivity  of  the  future  gravitational  wave  experiment  on  the  Cassini  mission.  Included  in  Table 
I  is  also  the  spectral  density  of  the  noise  of  a  crystal-driven  ultra-stable  oscillator  (USO). 

If  dual  radio  frequencies  in  the  uplink  and  downlink  are  used,  then  the  frequency  fluctuations 
due  to  the  interplanetary  plasma  can  be  entirely  removedl11!.  We  will  refer  to  this  configuration 
as  MODE  I.  If  only  one  frequency  is  adopted  instead,  which  we  will  assume  to  be  Ka-Band 
(32  GHz),  we  will  refer  to  this  configuration  as  MODE  II.  Ka-Band  is  planned  to  be  used  on 
most  of  the  forthcoming  NASA  missions,  and  will  be  implemented  on  the  ground  antennas  of 
the  Deep  Space  Network  (DSN)  by  the  year  1999  for  the  Cassini  mission. 

In  Figure  2  we  plot  the  r.m.s.  <?(/*)  of  the  noise  as  a  function  of  the  frequencies  /*  (k  = 
1,2,3, ....),  assuming  that  an  interplanetary  spacecraft  is  out  to  a  distance  L  =  1.0  AU.  For  this 
configuration  the  fundamental  frequency  of  the  xylophone  (Equation  (11))  is  equal  to  5.0  x  10-4 
Hz.  A  complete  analysis  covering  configurations  with  spacecraft  at  several  other  distances  is 
given  in!2!. 

The  MODE  I  configuration  is  represented  by  two  curves,  depending  on  whether  an  atomic  clock 
(circles)  or  a  USO  (squares)  is  operated  on  board  the  spacecraft.  Sensitivity  curves  for  the 
MODE  II  configuration  are  also  provided,  again  with  an  atomic  clock  on  board  (up-triangles) 
or  a  USO  (down-triangles).  The  best  sensitivity  is  achieved  in  the  MODE  I  configuration, 
regardless  of  whether  an  atomic  clock  or  a  USO  is  operated  on  board  the  spacecraft  (circles 
and  squares  are  over  imposed).  This  is  because  the  amplitude  of  the  noise  of  the  clock  on 
board  is  reduced  by  a  factor  7tA/L/c  =  4.7  x  10~*  at  the  xylophone  frequencies.  At  /  =  10~3 
Hz  the  corresponding  r.m.s.  noise  level  is  equal  to  4.7  x  10“18,  and  it  increases  to  a  value  of 
5.7  x  10-18  at  /  =  10~2  Hz.  As  far  as  the  MODE  II  configuration  is  concerned,  the  r.m.s. 
noise  level  is  equal  to  7.9  x  10"18  at  /  =  10-3  Hz,  while  it  decreases  to  6.3  x  10~1S  at  /  =  10"2 
Hz.  This  is  due  to  the  fact  that  the  one-sided  power  spectral  density  of  the  fractional  frequency 
fluctuations  due  to  the  interplanetary  plasma  decays  as  f~2/3. 

In  Figure  3  we  turn  to  redshift  experiments  with  a  spacecraft  out  to  5  AU,  and  we  assume  an 
observing  time  of  40  hours.  This  example  can  be  considered  as  representative  of  a  spacecraft 
orbiting  the  planet  Jupiter.  The  reduction  factor  of  the  buffeting  noise  B(t)  (see  Equation 
(10))  is  now  equal  to 


^=5.5 


x  10~2  , 


(14) 
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and  the  xylophone  frequencies  are  given  by  the  following  relation 

/fc  =  2l±_/;  fc  =  1,2,3, ....  (15) 

We  also  have  assumed  that  when  plasma  calibration  is  not  implemented,  the  frequency  fluctu¬ 
ations  due  to  interplanetary  scintillation  are  estimated  at  opposition.  This  of  course  does  not 
represent  a  general  situation  but  only  an  example. 

The  best  sensitivity  is  achieved  in  the  MODE  I  configuration  and  by  using  an  atomic  clock  on 
board.  At  the  Fourier  frequency  /  =  2.0  x  10-4  Hz  the  sensitivity  is  equal  to  3.0  x  10-17,  and 
it  increases  slowly  at  higher  frequencies.  In  the  MODE  II  configuration  and  with  an  atomic 
clock  on  board  the  sensitivity  degrades  by  about  a  factor  of  three  with  respect  to  the  previous 
configuration.  If  a  USO  is  used  instead,  then  MODE  I  and  II  are  totally  equivalent,  since 
the  USO  is  the  dominant  noise  source.  In  this  case  at  f  —  2.0  x  10-4  the  sensitivity  of  the 
xylophone  is  equal  to  1.4  x  10-14,  and  improves  as  2  as  the  frequency  increases. 

CONCLUSIONS 

We  have  discussed  a  method  for  significantly  increasing  the  sensitivity  of  two  Doppler  tracking 
experiments,  namely  searches  for  gravitational  waves  and  measurements  of  the  redshift  effect. 
Our  method  relies  on  a  properly  chosen  linear  combination  of  the  one-way  Doppler  data 
recorded  on  board  with  those  measured  on  the  ground.  It  allows  us  to  remove  entirely  the 
frequency  fluctuations  due  to  the  troposphere,  ionosphere,  and  antenna  mechanical,  and  for 
a  spacecraft  that  is  tracked  for  forty  days  out  to  1  AU  in  search  for  gravitational  waves,  it 
reduces  by  almost  four  orders  of  magnitude  the  noise  due  to  the  on-board  clock.  For  a  redshift 
experiment  instead,  with  a  spacecraft  out  to  5  AU,  our  technique  allows  us  to  reduce  by  about 
two  orders  of  magnitude  the  noise  of  the  antenna  and  buffeting  of  the  spacecraft. 

The  experimental  technique  presented  in  this  paper  can  be  extended  to  a  configuration  with 
two  spacecraft  tracking  each  other  through  a  microwave  or  a  laser  link.  Future  space-based 
laser  interferometric  detectors  of  gravitational  wavesl1*!,  for  instance,  could  implement  this 
technique  as  a  backup  option,  if  failure  of  some  of  their  components  would  make  the  normal 
interferometric  operation  impossible. 

As  a  final  note,  a  method  similar  to  the  one  presented  can  be  used  in  all  those  radio  science 
experiments  in  which  one-way  and  two-way  spacecraft  Doppler  measurements  are  used  as 
primary  data  set.  We  will  analyze  the  implications  of  the  sensitivity  improvements  that  this 
technique  will  provide  for  direct  measurements  of  the  following  quantities  such  as  searches  for 
possible  anisotropy  in  the  velocity  of  light,  measurements  of  the  Parameterized  Post-Newtonian 
parameters,  measurements  of  the  deflection  and  time  delay  by  the  sun  in  radio  signals,  and 
occultation  experiments.  This  research  is  in  phase  of  development,  and  will  be  the  subject  of 
a  forthcoming  paper. 
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Figure  1. 

Block  diagram  of  the  radio  hardware  at  the  ground  antenna  of  the  NASA  Deep  Space  Net¬ 
work  (DSN)  and  on  board  the  spacecraft  (S/C),  that  allows  the  acquisition  and  recording 
of  the  four  Doppler  data  F|(t),  Fj(/),  and  Sj(t). 
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Table  I 

List  of  the  noise  sources  entering  into  the  combined  Doppler  response  y(t).  The  .Allan 
deviation  at  a  given  integration  time  r  is  a  statistical  parameter  for  describing  frequency 
stability.  It  represents  the  root-mean-squared  expectation  value  of  the  random  process 
associated  with  the  fractional  frequency  changes,  between  time-contiguous  frequency  mea¬ 
surements,  each  made  over  time  intervals  of  duration  r.  The  numbers  provided  in  this 
table  are  taken  from  the  Riley  et  si.  report  [11]- 
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o  No  Plasma  &  Atomic  Clock 
□  No  Plasma  &  USO 
a  Plasma  &  Atomic  Clock 
v  Plasma  &  USO 


GRAVITATIONAL  WAVE  SEARCHES  (L  =  1  AU) 


Figure  2. 

The  r.m.s.  noise  level  as  a.  junction  of  the  frequencies  /t  (fc  =  1,2,3 . ),  estimated 

for  a  spacecraft  that  is  out  to  a  distance  L  =  1.0  AU  searching  for  gravitational  waves. 
Sensitivity  figures  for  the  four  distinct  configurations  are  represented  with  four  different 
symbols.  Circles  represent  r,m.s.  values  as  functions  of  the  xylophone  frequency  when 
plasma  frequency  fluctuations  are  totally  removed  and  an  atomic  standard  is  on  hoard. 
Squares  are  sensitivity  values  again  with  plasma  calibration,  but  now  with  a  USO  on  the 
spacecraft.  Up-triangles  are  used  for  representing  sensitivity  figures  affected  by  plasma 
noise  at  Ka-Band,  and  an  atomic  dock  is  used  on  board.  Down-triangles  assume  a  plasma 
noise  at  Ka-Band,  but  now  with  a  USO  on  board. 
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o  No  Plasma  &  Atomic  Clock 
□  No  Plasma  &  USO 
A  Plasma  &  Atomic  Clock 
v  Plasma  &  USO 


RED  SHIFT  EXPERIMENTS  (L  =  5  AU) 


Figure  3. 

The  r.m.s.  noise  level  as  a  function  of  the  frequencies  ft  (k  =  1,2, 3, estimated  for  a 
spacecraft  that  is  out  to  a  distance  L  =  5.0  AU  and  tracked  during  a  red  shift  experiment. 


